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ABSTRACT 

This paper presents new deep and wide narrow-band surveys undertaken with UKIRT, Subaru 
and the VLT; a unique combined effort to select large, robust samples of Ha star-forming 
galaxies at z = 0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 
9.2 and 10.6 Gyrs) in a uniform manner over ~ 2 deg^ in the COSMOS and UDS fields. 
The deep multi-epoch Ha surveys reach a matched 3 a flux limit of ^ 3M0yr~^ out to 
z = 2.2 for the first time, while the wide area and the coverage over two independent fields 
allow to greatly overcome cosmic variance and also assemble large samples of more luminous 
galaxies. Catalogues are presented for a total of 1742, 637, 515 and 556 Ha emitters, robustly 
selected (E>3, EWo(Ha)>25 A) at 2; = 0.40, 0.84, 1.47 and 2.23, respectively, and used to 
determine the Ha luminosity function at those epochs and its evolution. The faint-end slope 
of the Ha luminosity function is found to be a = — 1.60±0.08 over z = — 2.23, showing no 
significant evolution. The characteristic luminosity of SF galaxies, L^^, evolves significantly 
as log L^ai^) = OAbz + log L^^q. This is the first time Ha has been used to trace SF 
activity with a single homogeneous survey dXz = 0.4 — 2.23. Overall, the evolution seen in the 
Ha luminosity function is in good agreement with the evolution seen using inhomogeneous 
compilations of other tracers of star formation, such as FIR and UV, jointly pointing towards 
the bulk of the evolution in the last 1 1 Gyrs being driven by a similar star-forming population 
across cosmic time, but with a strong luminosity increase from z~0toz~2.2. Our uniform 
analysis allows to derive the Ha star formation history of the Universe, showing a clear rise 
up to z 2.2, for which the simple parametrisation logpsFR = — 2.1/(z + 1) is a good 
approximation for z < 2.2. The results reveal that both the shape and normalisation of the Ha 
star formation history are consistent with the measurements of the stellar mass density growth, 
confirming that our Ha cosmic star formation history is tracing the bulk of the formation of 
stars in the Universe for 2: < 2.23. The star formation activity over the last ^11 Gyrs is 
responsible for producing ~ 95% of the total stellar mass density observed locally today, with 
about half of that being assembled at 1.2 < z < 2.2, and the other half at z < 1.2. 

Key words: galaxies: high-redshift, galaxies: luminosity function, cosmology: observations, 
galaxies: evolution. 



This work is based on observations obtained using the Wide Field CAM- 
era (WFCAM) on the 3.8m United Kingdom Infrared Telescope (UKIRT), 
as part of the High-redshift(Z) Emission Line Survey (HiZELS; U/CMP/3 
and U/lOB/07). It also relies on observations conducted with HAWK-I on 
the ESO Very Large Telescope (VLT), program 086.7878.A, and observa- 
tions obtained with Suprime-Cam on the Subaru telescope (S10B-144S). 
t NOVA Fellow; E-mail: sobral@strw.leidenuniv.nl 



1 INTRODUCTION 

Star formation activity in galaxies has been decreasing signifi- 
cantly over time. Such observational evidence can now be seen very 
clearly by studies looking at how the specific star formation rate 
( |Salim et al.|2007l [Gonzalez et al.|2010l|Dutton et al.|2010| ), or the 
star formation rate density (psfr) vary as a function of look-back 
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time. Nevertheless, while surveys reveal that psfr rises steeply out 
to at least z ^ I (e.g. [Lilly et al.|1996[p3pkins & Beacom|2006| ), 
determining the redshift where psFR might have peaked at z > 1 
is still an open problem. This is because the use of different tech- 
niques/indicators (affected by different biases, dust extinctions and 
with different sensitivities - and that can only be used over lim- 
ited redshift windows) results in a very blurred and scattered under- 
standing of the star formation history of the Universe. Other prob- 
lems/limitations result from the difficulty of obtaining both large- 
area, large-sample, clean and deep observations (to overcome both 
cosmic variance, and avoid large extrapolations down to faint lumi- 
nosities). 

One route towards achieving significant progress is through 
the use of the narrow-band technique to undertake large, sensitive 
surveys and trace emission lines from the optical into the near-IR. 
These are particularly powerful, as they can take advantage of large 
field-of-view detectors with narrow-band filters and select sources 
purely on the basis of their line emission. While there is a wide 
range of emission lines which can trace star formation. Ha is by 
far the best candidate to do so over 2; < ^ as it provides a sensi- 
tive census of star formation activity, is very well-calibrated, and is 
typically affected by only 1 mag of extinction at 2; ^ — 1.5 (e.g. 
[Gilbank et al.||2010| [Sobral et al.||2012| ), contrary to other lower 
wavelength emission lines. Longer wavelength lines such as the 
Paschen series lines are less affected by dust extinction, but they 
are intrinsically fainter than Ha (e.g. Paa is intrinsically ^ 10 x 
weaker than Ha for a typical star-forming galaxy) and much harder 
to observe, even for moderate redshifts. 

Ha surveys have been carried out by many authors (e.g. 
[Bunker et al.|1995|[Malkan et al.|1995 l, but they initially resulted 
in a relatively low number of sources for z > 0.5 surveys. Fortu- 
nately, the development of wide field near-IR detectors has recently 
allowed a significant increase in success: at 2; ^ 2, narrow-band 
surveys such as [Moorwood et al.] { |2000 ), which could only detect 
a handful of emitters, have been rapidly extended by others, such 
as [Geach et al.[ ( [2008| ), increasing the sample size by more than 
an order of magnitude. Substantial advances have also been ob- 
tained at z - 1 (e.g. [Villar et al.|2008[ [Sobral et al.|2009|[Ly et al.[ 
[201 1). Other Ha surveys have used dispersion prisms on HST to 
make progress (e.g. McCarthy et al.[1999|| Yan et al.[199"9l [Hopkins] 
[et al.[200^ , and there is promising work being conducted using the 
upgraded WFC3-grism (e.g. WISP or 3D-HS T; [Atek et al.[[20T0l 
[Straughn et al.|2011[[van Dokkum et al.|2011| ). 

HiZELS, the High-redshift(Z) Emission Line Surve}[^('Geach' 
[et al.|2008[ [Sobral et al.|2009 l|2UT2l hereafter S09 and SI 2) is a 
Campaign Project using the Wide Field CAMera (WFCAM) on 
the United Kingdom Infra-Red Telescope (UKIRT) and exploits 
specially-designed narrow-band filters in the J and H bands (NBj 
and NBh), along with the H2S(1) filter in the K band (hereafter 
NBk), to undertake panoramic, moderate depth surveys for line 
emitters. HiZELS is primarily targeting the Ha emission line red- 
shifted into the near-infrared at z = 0.84, z = 1.47 and z = 2.23 
(see [Best et al.|2010| ), while the NBj and NBh filters also provide 

^ There may be possible to extend such studies to even higher redshifts: 
[Shim et all] ( [201 1) provided evidence that Spitzer fluxes can be highly con- 
taminated by Ha emission at even higher redshifts (z ~ 4), and thus such 
excess - if robustly measured, can in principle be used to identify Ho; emit- 
ters up to even earlier epochs in the Universe, jwst/nircam will certainly 
open a new window to explore Ho; emission at such high redshifts. 
^ For more details on the survey, progress and data releases, see 
http://www.roe.ac .uk/ifa/HiZELS/| 



Table 1. Narrow-band filters used to conduct the multi-epoch surveys for 
Hq; emitters, indicating the central wavelength (|Ltm), full width at half max- 
imum (FHWM), the redshift range for which the Ho; line is detected over 
the filter FWHM, and the corresponding volume (per square degree) sur- 
veyed (for the Hq; line). Note that the NB921 filter provides an [On] survey 
which precisely matches the Ho; z = 1.47 survey, and also a [Oiii] survey 
which broadly matches the 2; = 0.84 Ha survey. The NBj and NBh filters 
also provide [On] 3727 and [Oiii] 5007 surveys, respectively, which match 
the z = 2.23 NBk Ha survey 



NB filter 


Ac 


FWHM 


zHa 


Volume (Hq;) 




(M-m) 


(A) 




(10^ Mpc^ deg-2) 


NB921 


0.9196 


132 


0.401±0.010 


5.13 


NBj 


1.211 


150 


0.845±0.015 


14.65 


NBh 


1.617 


211 


1.466±0.016 


33.96 


NBk 


2.121 


210 


2.231±0.016 


38.31 


HAWK-I H2 


2.125 


300 


2.237±0.023 


5.47 



[On] 3727 and [Olll] 5007 surveys slI z ^ 2.23, matching the NBk 
Ha coverage at the same redshift. 

One of the main aims of HiZELS is to provide measurements 
of the evolution of the Ha luminosity function from z — 0.0 to 
z = 2.23 (but also other properties, such as clustering, environ- 
ment and mass dependences; c.f. [Sobral et al . 2010 2011). The 
first results (Geach et al. 2008 S09; S12) indicate that the Ha lu- 
minosity function evolves significantly, mostly due to an increase 
of about one order of magnitude in Lh^, the characteristic Ha lu- 
minosity, from the local Universe to 2; = 2.23 (S09). In addition, 
Sobral et al. ( 201 1) found that at z = 0.84 the faint-end slope of the 
luminosity function (a) is strongly dependent on the environment, 
with the Ha luminosity function being much steeper in low density 
regions and much shallower in the group/cluster environments. 

However, even though the progress has been quite remarkable, 
significant issues remain to be robustly addressed for a variety of 
reasons. For example, is the faint-end slope of the Ha luminosity 
function (a) becoming steeper from low to high redshift? Results 
from [Hayes et"ar] { [2010| ) point towards a steep faint-end slope at 
z > 2. However, Hayes et al. did not sample the bright end, and 
have only targeted one single field over a relatively small area, and 
thus cosmic variance could play a huge role. [Tadaki et al.[ { [201 1| ) 
find a much shallower a at 2; ^ 2 using Subaru. Furthermore, mea- 
surements so far rely on different data, obtained down to different 
depths and using different selection criteria. Additionally, differ- 
ent ways of correcting for completeness (c.f. for example Ly et al.] 
[201 1| ), filter profiles or contamination by the [Nii]aa6548,6583.6 
lines can also lead to significant differences. How much of the evo- 
lution is in fact real, and how much is a result of different ways 
of estimating the Ha luminosity function? This can only be fully 
quantified with a completely self-consistent multi-epoch selection 
and analysis. Another issue which still hampers the progress is 
overcoming cosmic variance and probing a very wide range of en- 
vironments and stellar masses at z > 1. Large samples of homo- 
geneously selected star-forming galaxies at different epochs up to 
z > 2 would certainly be ideal to provide strong tests on our un- 
derstanding of how galaxies form and how they evolve. 

In order to clearly address the current shortcomings and pro- 
vide the data that is required, we have undertaken by far the largest 
area, deep multi-epoch narrow-band Ha surveys over two differ- 
ent fields. By doing so, both faint and bright populations of equally 
selected Ha emitters at z = 0.4, 0.84, 1.47 and 2.23 have been 
obtained, using 4 narrow band filters (see Figure 1 and Table 1). 
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Table 2. Observation log of all the narrow-band observations obtained over the COSMOS and UDS fields, taken using WFCAM on UKIRT, HAWK-I on the 
VLT, and Suprime-Cam on Subaru, during 2006-201 1. Limiting magnitudes (3 a, all in Vega) are the average of each field/pointing, based on the measurements 
of 10^-10^ randomly placed 2" apertures in each frame. 



Field 


Band 


R.A. 


Dec. 


Int. time 


FHWM 


Dates 


miim(Vega) 




(filter) 


(J2000) 


(J2000) 


(ks) 


n 




(3C7) 


COSMOS- 1 


NB921 


09 59 23 


+02 3029 


2.9 


0.9 


2010 Dec 9 


24.4 


COSMOS-2 


NB921 


100134 


+02 3029 


2.9 


0.9 


2010 Dec 9 


24.5 


COSMOS-3 


NB921 


09 59 23 


+02 0416 


2.9 


0.9 


2010 Dec 9 


24.5 


rOSMOS-4 


NB921 


1001 34 


^02 04 1 6 


2.9 


0.9 


201 Dec 9 


24.5 


UKIDSS-UDS C 


NB921 


0218 00 


-05 0000 


30.0 


0.8 


2005 Oct 29, Nov 1, 2007 Oct 11-12 


26.6 


UKIDSS-UDS N 


NB921 


0218 00 


-04 35 00 


37.8 


0.9 


2005 Oct 30,31, Nov 1, 2006 Nov 18, 2007 Oct 11,12 


26.7 


UKIDSS-UDS S 


NB921 


0218 00 


-05 25 00 


37.1 


0.8 


2005 Aug 29, Oct 29, 2006 Nov 18, 2007 Oct 12 


26.6 


UKIDSS-UDS E 


NB921 


021947 


-05 0000 


29.3 


0.8 


2005 Oct 31, Nov 1, 2006 Nov 18, 2007 Oct 11,12 


26.6 


UKIDSS-UDS W 


NB921 


021613 


-05 0000 


28.1 


0.8 


2006 Nov 18, 2007 Oct 11,12 


26.0 


COSMOS-NW(l) 


NBj 


100000 


-h021030 


19.7 


0.8 


2007 Jan 14-16 


22.0 


COSMOS-NE(2) 


NBj 


100052 


-F021030 


23.8 


0.9 


2006 Nov 10; 2007 Jan 13-14 


22.0 


COSMOS-SW(3) 


NBj 


100000 


-h0223 44 


18.9 


0.9 


2007 Jan 15-17 


22.0 


COSMOS-SE(4) 


NBj 


100053 




17.1 


1.0 


2007 Jan 15, 17; Feb 13, 14, 16 


21.9 


UKIDSS-UDS NE 


NBj 


0218 29 


-04 5220 


20.9 


0.8 


2007 Oct 21-23 


22.0 


UKIDSS-UDS NW 


NBj 


0217 36 


-04 5220 


22.4 


0.9 


2007 Oct 20-21 


22.1 


UKIDSS-UDS SE 


NBj 


0218 29 


-05 05 53 


19.6 


0.9 


2007 Oct 23, 24 


22.0 


UKIDSS-UDS SW 


NBj 


0217 38 


-05 05 34 


22.4 


0.8 


2007 Oct 19, 21 


22.0 


V_.wOiVlvJv3-i> VV i^i^ 






-1-09 1 n 


iZ. J 


1 n 


9000 Inp'K 97- Mnr 1 9 
Z.\j\jy -TCD z/, iviar i-Z 


91 1 
Zi. 1 


COSMOS-NE(2) DEEP 


NBh 


100052 


+021030 


107.0 


0.9 


2009 Feb 28; Apr 19; May 22; 201 1 Jan 26-30 


22.2 


COSMOS-SW(3) 


NBh 


100000 


+02 23 44 


14.0 


0.7 


2010 Apr 2 


21.0 


COSMOS-SE(4) 


NBh 


100053 


+02 23 44 


18.1 


1.0 


2009 Mar 2; Apr 30; May 22; 2010 Apr 3 


20.8 


COSMOS-A 


NBh 


100001 


+02 3653 


12.6 


1.0 


2010 Apr 9; 2011 Jan 25 


21.1 


COSMOS-B 


NBh 


100054 


+02 36 30 


12.6 


0.9 


2010 Apr 8-9 


21.0 


COSMOS-C 


NBh 


100001 


+015710 


13.0 


0.8 


2010 Apr 6-8 


20.8 


COSMOS-D 


NBh 


100052 


+015715 


14.2 


0.8 


2010 Apr 7-8 


21.0 


COSMOS-E 


NBh 


09 59 07 


+02 23 44 


12.6 


0.8 


2010 Apr 3-6 


20.9 


COSMOS-F 


NBh 


09 59 07 


+021030 


12.6 


0.7 


2010 Apr 4 


21.1 


COSMOS-G 


NBh 


100146 


+02 23 44 


14.5 


0.8 


2010 Apr4&9 


21.0 


COSMOS-H 


NBh 


100148 


+021051 


14.0 


0.8 


2010 Apr 6 & 9 


20.8 


UKIDSS-UDS NE 


NBh 


0218 29 


-04 5220 


18.2 


0.9 


2008 Sep 28-29; 2009 Aug 16-17; 2010 Jul 22 


21.3 


UKIDSS-UDS NW 


NBh 


0217 36 


-04 5220 


18.0 


0.9 


2008 Sep 25, 29; 2010 Jul 18,22 


20.9 


UKIDSS-UDS SE 


NBh 


0218 29 


-05 05 53 


25.2 


0.8 


2008 Sep 25, 28-29; 2009 Aug 16-17 


21.5 


TIKIDSS-TIDS SW 


NRu 


02 17 38 


—05 05 34 


19.6 


0.9 


200R Oct-Nov 2009 An? 16-17- 2010 Tnl 23 


21.3 


COSMOS-NW(l) 


NBk 


100000 


+021030 


24.3 


0.9 


2008 May 11; 2009 Feb 27 


21.0 


COSMOS-NE(2) DEEP 


i"N13K 


1 00 52 


+02 1 30 


62.5 


1.0 


2006 Mav 20-26* Dec 5 1 4-20- 2008 Mar 6-9 


21.3 


COSMOS-SW(3) 


NBk 


100000 


+02 23 44 


16.5 


0.9 


2006 May 20-26, Dec 5, 14-20; 2009 May 20 


20.6 


COSMOS-SE(4) 


NBk 


100053 


+02 23 44 


18.0 


0.9 


2006 Nov 13-17 


20.7 


COSMOS-A 


NBk 


100001 


+02 3653 


20.0 


1.0 


2011 Mar 19-26 


20.9 


COSMOS-B 


NBk 


100054 


+02 36 30 


20.0 


0.9 


2011 Mar 25-26 


20.7 


COSMOS-C 


NBk 


100001 


+015710 


26.7 


0.8 


2011 Mar 27, 30; Apr 3-5, 16-18 


20.9 


COSMOS-D 


NBk 


100052 


+015715 


6.7 


0.8 


2011 Mar 30; Jun 6 


20.2 


COSMOS-E 


NBk 


09 59 07 


+02 23 44 


20.0 


0.8 


2011 Mar 30; May 18-23, 30 


20.7 


UKIDSS-UDS NE 


NBk 


0218 29 


-04 5220 


19.2 


0.8 


2005 Oct 18; 2006 Nov 13-14 


20.7 


UKIDSS-UDS NW 


NBk 


0217 36 


-04 5220 


20.0 


0.9 


2006 Nov 11 


20.8 


UKIDSS-UDS SE 


NBk 


0218 29 


-05 05 53 


18.7 


0.8 


2006 Nov 15-16 


20.7 


UKIDSS-UDS SW 


NBk 


0217 38 


-05 05 34 


23.6 


0.8 


2007 Sep 30 


20.9 


COSMOS-HAWK-I 


H2 


100000 


+021030 


17.1 


0.9 


2009 Apr 10, 14-15, 18; May 13-14 


21.5 


UKIDSS-UDS-HAWK-I 


H2 


0217 36 


-04 5220 


17.1 


1.0 


2009 Aug 16, 19, 24, 27 


21.7 



This paper presents the narrow-band imaging and results obtained 
with the NB921, NBj, NBh, NBk and H2 filters on Subaru, the 
United Kingdom InfraRed Telescope (UKIRT) and the Very Large 
Telescope (VLT), over a total of ~ 2 deg^ in the Cosmological 
Evolution Survey (COSMOS; [Scoville et all2007^ ) and the SXDF 
Subaru-XMM-UKIDSS Ultra Deep Survey (UDS; |LawFence et al.| 
[2007] ) fields. 



The paper is organised as follows: §2 describes the observa- 
tions, data reduction, source extraction, catalogue production, se- 
lection of line emitters, and the samples of Ha emitters. In §3, after 
estimating and applying the necessary corrections, the Ha luminos- 
ity functions are derived at z = 0.4, 0.84, 1.47 and 2.23, together 
with an accurate measurement of their evolution at both bright and 
faint ends. In §4 the star formation rate density at each epoch is 
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Figure 1. The broad- and narrow-band filter profiles used for the analysis. 
The narrow-band filters in the z' , J, H and K bands (typical FWHM of 
^ 100 - 200 A) trace the redshifted Ha line at z = 0.4, 0.84, 1.47, 2.23 
very effectively, while the (scaled) broad-band imaging is used to estimate 
and remove the contribution from the continuum. Note that because the fil- 
ters are not necessary located at the center of the respective broad-band 
transmission profile, very red/blue sources can produce narrow-band ex- 
cesses which mimic emission lines; that is corrected by estimating the con- 
tinuum colour of each source and correcting for it. 



also evaluated and the star formation history of the Universe is pre- 
sented. §4 also discusses the results in the context of galaxy forma- 
tion and evolution in the last 1 1 Gyrs, including the inferred stel- 
lar mass density growth. Finally, §5 presents the conclusions. An 
Ho = 70km s~^ Mpc~^, Q.m — 0.3 and — 0.7 cosmology is 
used. Narrow-band magnitudes in the near-infrared and the associ- 
ated broad-band magnitudes are in the Vega system, except when 
noted otherwise (e.g. for colour-colour selections). NB921 and z' 
magnitudes are given in the AB system. 



2 DATA AND SAMPLES 

2.1 Optical NB921 imaging with Subaru 

Optical imaging data were obtained with Suprime-Cam using the 
NB921 narrow-band filter. The NB921 filter is centered at 9196 A 
with a FWHM of 132 A. The COSMOS field was observed in ser- 
vice mode in December 2010 with four different pointings covering 
the central 1.1 deg^. Total exposure times were 2.9 ks per pointing, 
composed of individual exposures of 360 s dithered over 8 different 
positions. Observations are detailed in Table 2. The UDS field has 
also been observed with the NB921 filter (see |Ouchi et al. 12010] ), 
and these data have been extracted from the archive. Full details of 
the data reduction and catalogue production of the UDS data were 
presented by S12 and the same approach was adopted for the COS- 
MOS data. In brief, all the raw NB921 data were reduced with the 
Suprime-Cam Deep field REDuction package (SDFRED, |Yagi et al"] 
|2002[ |Ouchi et aL]|2004| ) and IRAF. The combined images were 
aligned to the public z^-band images of Subaru-XMM Deep Sur- 
vey or the COSMOS field and PSF matched (FWHM= 0.9'0- The 
NB921 zero points were determined using z' data, so that the {z'- 
NB921) colours are consistent with a median of zero for z' between 
19 and 21.5- where both NB921 and z' images are unsaturated and 
have very high signal-to-noise ratios. 
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Figure 2. The survey strategy used to cover the COSMOS field. The cen- 
tral pointings (1,2,3,4; see Table 2) were complemented with further point- 
ings (A to H) to both increase the surveyed area, and increase the exposure 
time in the central area. The region delimited by the dashed line shows 
the NB921 coverage obtained in COSMOS. See Table 2 for details on the 
pointings, including exposure times. 



Source detection and photometry were performed using S EX- 
TRACTOR ( [Bertin & Arnouts 11996] ). Sources were detected on each 
individual NB921 image and magnitudes measured with 2" and 3'^ 
diameter apertures. The 3^^ apertures are used to select and mea- 
sure Ha line fluxes: at z = 0.4 the 3^^ apertures measure the same 
physical area as 1" apertures dXz — 0.8, 1.47, 2.23 16kpc), as- 
suring full consistency. The 2" apertures are used to measure emis- 
sion lines from sources at higher redshift ([Oil] z — 1.47, to 
match the NBh Ha measurement at the same redshift, and [Olll] 
at z = 0.84 to match the NBj Ha survey). The average NB921 3 a 
limiting magnitudes (in 2" apertures) are given in Table[2] 



2.2 Near-infrared imaging with UKIRT 

The COSMOS and UKIDSS UDS fields were observed with WF- 
CAM on UKIRT as summarised in Table 2, using the NBj, NBr 
and NBk narrow-band filters, with central wavelengths and FWHM 
given in Table 1. Observations were conducted over 2006-2011, 
covering 1.6 deg^ (NBj), 2.34 deg^ (NBh) and 2.0 deg^ (NBk) 
over the COSMOS and the UDS fields (see Table 2). 

For the COSMOS field, an initial 0.8 deg^ coverage obtained 
with 4 WFCAM pointings was complemented in NBh (and NBk) 
by 8 (5) further WFCAM pointings, macro-jittered to obtain a com- 
bined 1.6 (1.3) deg^ coverage with increasing exposure time per 
pixel towards the centre of the field (see Figure [2]). Part of the cen- 
tral region (^ 0.2 deg^) benefits further from some significant ex- 
tra deep data (see Table 2), leading to a much higher total exposure 
time. 

A dedicated pipeline has been developed for HiZELS 
(PfHiZELS, c.f. S09 for more details). The pipeline has been mod- 
ified and updated since S09 mostly to 1) improve the flat-fieldin£] 

^ The improvements in the fiat-fielding are obtained by stacking all second- 
pass flattened frames per field and producing source masks on the stacked 
images. Masks are then used to produce third pass fiats using all the frames 
in the jitter sequence except the frame being flattened. The third-pass flat- 
tened frames are then stacked again, and the procedure is repeated another 
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and 2) provide more accurate astrometric solutions for each individ- 
ual frame which result in a more accurate stacking The updated 
version of the pipeline (PfHiZELS2012) has been used to reduce all 
UKIRT narrow band data (NBj, NBh and NBk), including those 
already presented in previous papers. This approach guarantees a 
complete self-consistency and takes advantage of the improved re- 
duction which, in some cases, is able to go deeper by ^ 0.2 mag 
when compared to the data reduced by the previous version of the 
pipeline. 

For the COSMOS field, due to the need to co-add frames taken 
with different WFCAM cameras (due to the survey strategy, see 
Figure |2]), IRAF is used to distort-correct all final frames (in addi- 
tion to individual frames being corrected prior to combining) before 
co-adding different fields; the typical rms is 0.1 — 0.2^^ However, at 
the highest radial distances (r > 1000 pix) from the centre of the 
images, even small residual distortions (sometimes with opposite 
signs due to the need of combining different WFCAM cameras) 
can lead to an elongated PSF and blurring of sources. The com- 
bined frames were analysed using IRAF and it was found that for 
r < 900 (pix, from the centre of the images obtained with different 
WFCAMs) the PSF/ellipticity remained unchanged due to stack- 
ing, but that for higher radial distances the PSF became increas- 
ingly elongated, particularly in the comers. In order to avoid intro- 
ducing errors in the photometry and source position, regions with 
overlapping WFCAM cameras are treated in the following way: i) 
for r < 900 pix, the stacked images are used and ii) for r > 900 
pix the deepest individual image is used. 

Narrow-band images were photometrically calibrated (inde- 
pendently) by matching 100 stars per frame with J, H and K 
between the 12th an 16th magnitudes from the 2MASS All-Sky 
catalogue of Point Sources (Cutri et al. 2003 ) which are unsat- 
urated in the narrow-band images. WFCAM images are affected 
by cross-talk and other artifacts caused by bright stars: accurate 
masks are produced in order to reject such regions. Sources were 
extracted using S Extractor ( |Bertin & Arnouts|[l996| ), making 
use of the masks. Photometry was measured in apertures of 2" di- 
ameter which at z = 0.8 — 2.2 recover Ha fluxes over ^ 16kpc. 
The average 3 a depths of the entire set of NB frames vary signifi- 
cantly, and are summarised in Table 2. The total numbers of sources 
detected with each filter are given in Table 3. 

2.3 Near-infrared H2 imaging with HAWK-I 

The UKIDSS UDS a nd COSMOS fields were ob served with the 
HAWK-I instrument ( |Pirard et al.|2004||CasaH et al.„2006 ) on the 
VLT during 2009. A single dithered pointing was obtained in each 
of the fields using the H2 filter, characterised by Ac = 2.124 |^m 
and SX = 0.030 ^im (note that the filter is slightly wider than that 
on WFCAM). Individual exposures were of 60 s, and the total ex- 
posure time per field is 5 hours. Table 2 presents the details of the 
observations and depth reached. 

Data were reduced using the HAWK-I ESO pipeline recipes, 
by following an identical reduction scheme/procedure to the WF- 
CAM data. The data have also been distortion corrected and as- 
trometrically calibrated before combining, using the appropriate 

time. This procedure is able to both mask many sources which are unde- 
tected in individual frames out of the flats, but particularly to mask bright 
sources much more effectively, as the stacking of all images reveals a wider 
distribution of flux from those sources. 

IRAF is used to distort correct the frames and obtain the best astrometric 
solution for each frame, always assuring that the flux is conserved. 



pipeline recipes. After combining all the individual reduced frames 
it is possible to obtain a contiguous image of ^ 7.5 x 7.5 arc 
min^ in each of the fields. There are, nonetheless, small regions 
with slightly lower exposure time per pixel in regions related with 
chip gaps at certain positions. Because of the availability of the 
very wide WFCAM imaging, regions in the HAWK-I combined 
images for which the exposure time per pixel is < 80% of the to- 
tal are not considered. Frames are photometrically calibrated using 
2MASS as a first pass, and then using UDS and COSMOS Ks cal- 
ibrated images to guarantee a median colour (K-NBr) for all 
magnitudes probed, as this procedure provides a larger number of 
sources. Similarly to the procedure used for WFCAM data, sources 
were extracted using S Extractor and photometry was measured 
in apertures of 2" diameter. 

2.4 Narrowband excess selection 

In order to select potential line emitters, broad-band (BB) imag- 
ing is used in the z\ J, H and Ks bands to match narrow-band 
(NB) imaging in the NB921, NBj, NBh and NBk/H2, respec- 
tively. Count levels on the broad-band images are scaled down to 
match the counts (of 2MASS sources) for each respective narrow- 
band image, in order to guarantee a median zero colour, and a com- 
mon counts-to-magnitude zero point. Sources are extracted from 
BB images using the same aperture sizes used for NB images 
and matched to the NB catalogue with a search radius of < 0.9^^ 
Note, however, that none of the narrow-band filters fall at the centre 
of the broad-band filters (see Figure 1). Thus, objects with signifi- 
cant continuum colours will not have BB — NB = 0; this can be 
corrected with broad-band colours (c.f. SI 2), in order to guarantee 
that BB — NB distribution is centred on and has no dependence 
on continuum broad-band colours. Average colour correction^ are 
given by: 

(z^-NB921)ab = (z^-NB921)o,ab-0.05(Jab-z^b)-0.04 
(J - NBj) = (J - NBj)o - 0.09(z:vB - Jab) + 0.11 
(H - NBh) = (H - NBh)o + 0.07(Jab - Hab) + 0.06 
(K - NBk) = (K - NBk)o - 0.02(Hab - Kab) + 0.04 

Potential line emitters are then selected according to the sig- 
nificance of their (BB — NB) colour, as they will have (BB — 
NB) > 0. True emitters are distinguished from those with positive 
colours due to the scatter in the magnitude measurements by quan- 
tifying the significance of the narrowband excess. The parameter S 
(see B unker et al.,1995) quantifies the excess compared to the ran- 
dom scatter expected for a source with zero colour, as a function of 
narrow-band magnitude (see e.g. S09), and is given by: 

I _ IQ-0A{BB-NB) 

where ZP is the zero point of the NB (and BB, as those have been 
scaled to have the same ZP as NB images), rap is the aperture 
radius (in pixels) used, and ctnb and ctbb are the rms (per pixel) of 
the NB and BB images, respectively. 

Here, potential line emitters are selected if S > 3.0. The 
spread on the brighter end (narrow-band magnitudes which are not 
saturated, but for which the scatter is not affected by errors in the 

^ Sources for which one of the broad-band colours is not available are as- 
signed the median correction. The median corrections are: —0.04, +0.07, 
+0.05 and +0.03 for NB921, NBj, NBh and NBk respectively 
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Figure 3. Photometric redshift distributions (peak of the probabiHty distribution function for each source) for the NB921 Hne emitter candidates (upper left). 
NB J (upper right), NBh (lower left) and NBk (lower right). All distributions peak at redshifts which correspond to strong emission lines, corresponding to 
Wot, H/3, [Olll] and [On] emitters at various redshifts depending on the central wavelength of each narrow-band filter. Other populations of emitters are also 
found, such as Paschen-lines. The dotted lines indicate the redshift for which an emission line is detectable by the narrow-band filters. Lya emitter candidates 
in the NB921 dataset (e.g. 18 sources in COSMOS) with photo-z of 2; ~ 6 — 7 are not shown. 



magnitude, i.e., much brighter than the limit of the images) is quan- 
tified for each data- set and frame, and the minimum {BB — NB) 
colour limit over bright magnitudes is set by the 3 x the standard 
deviation of the excess colour over such magnitudes (sb). A com- 
mon rest-frame EW limit of EWo = 25 A is applied, guaranteeing a 
limit higher than the 3 x dispersion over bright magnitudes in all 
bands. The combined selection criteria guarantees a clean selection 
of line emitters and, most importantly, it ensures that the samples of 
Ha emitters are selected down to the same rest-frame EW, allowing 
one to quantify the evolution across cosmic time. 



As a further check on the selection criteria, the original imag- 
ing data are used to produce BB and NB postage stamp images 
of all the sources. The BB is subtracted from the NB image leav- 
ing the residual flux. From visual inspection these residual images 
contain obvious narrow band sources and it is found that the re- 
maining flux correlates well with the catalogue significance, with 
significance being proportional to residual flux (/) as/^°-^^=^° °^^. 



2.5 The samples of NB line emitters 

Narrow-band detections below the estimated 3 a detection thresh- 
old were not considered. By using colour-colour diagnostics (see 
S 12), potential stars are identified in the sample and rejected as well 
(the small fraction varies from band to band; see Table 3). The sam- 
ple of remaining potential emitters (S > 3 & EWo > 25 A; see Ta- 
ble 3 for numbers) is visually checked to identify spurious sources, 
artifacts which might not have been masked, or sources being iden- 
tified in very noisy regions. Sources classed as spurious/artifacts are 
removed from the sample of potential emitters. The final samples 
of line emitters are then derived. 

As a further test of the reliability of the line emitter samples, it 
can be noted that since the HAWK-I observations are both deeper 
and obtained over a larger redshift slice (due to a wider filter pro- 
file) when compared to WFCAM, they should be able to confirm 
all NBk emitters over the matched area. This is confirmed, as all 
10 emitters which are detected with WFCAM in the matched area 
are recovered by HAWK-I data as well. 

The catalogues, containing all narrow-band emitter candi- 
dates, are presented in Appendix A. The catalogues provide IDs, 
coordinates, narrow-band and broad-band magnitudes, estimated 
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fluxes and observed EWs. Further details and information are avail- 
able on the HiZELS website. 

The photometric redshift (Photo-zs from [Ilbert e t al. 2009! 
[Cirasuolo et al.| |2010 ) distributions of the sources selected 
with the 4 narrow-band filters are presented in Figure [3] The 
photometric redshifts show clear peaks associated with Ha, 
H/3/[Oiii]aa4959,5007, and [Oii]a3727 (see Figure[3|, together with 
further emission lines such as Paschen lines and hya. Spectro- 
scopic redshift are also available for a fraction of the selected line- 
emitters^ 



Geach et al ||2008[ [van Breukelen et al.||2007[ |Ouchi et al.|[2008 



Small et al.|2008l|Ono et aL|2010|' 



2.6 Selecting Ha emitters 

Samples of Ha emitters at the various redshifts are selected using a 
combination of broad-band colours (colour-colour selections) and 
photometric redshifts (when available). Colour-colour separation of 
emitters are different for each redshift, and for some redshifts two 
sets of colour-colour separations are used to reduce contamination 
to a minimum. Additionally, spectroscopically confirmed sources 
are included, and sources confirmed to be other emission lines re- 
moved from the samples - but the reader should note that for all 
four redshifts the number of z-included sources missed by the se- 
lection is typically < 10, while the z-rejected sources are typically 
< 10 as well. Therefore, the decrease of available spectroscopic 
redshifts with redshift does not introduce any bias. 

Additionally, sources found to be line emitters in two (or three, 
for Ha emitters at z = 2.23) bands, making them strong Ha can- 
didates are also included in the samples, even if they have been 
missed by the colour-colour and photometric selection (although it 
is found that only very few real Ha sources are missed by the se- 
lection criteria). Table [Sjprovides the number of sources, including 
spectroscopically confirmed ones for each field at each redshift. 

2.6.1 Ha emitters at z = 0.4 

The selection of Ha emitters at z = 0.4 is done using the BRiK 
(B — R \s i — K) colour-colour selection (see Figure|4]and S12). 
In addition, all sources for which 0.35 < ^phot < 0.45 are also 
included in the sample. Spectroscopic redshifts (from zCOSMOS; 
[Lilly et al.|2009| ) are available for 38 sources. Thirty six sources are 
confirmed to be at z = 0.391 — 0.412, while 2 sources are [Nil] 
emitters. This implies a very high completeness of the sample, and 
a contamination of ^ 5 per cent over the entire sample. Contami- 
nants have been removed, and spectroscopic sources added. A total 
of 1742 Ha emitters at 2; = 0.4 are selected. 



2.6.2 Ha emitters at z — 0.84 

Sources are selected to be Ha emitters at 2; = 0.84 if 0.75 < 
^^phot < 0.95 or if they satisfy the BRiK (see Figure |4) S09) 
colour-colour selection fox z ^ 0.8 sources. Additionally, sources 
with 1.3 < Zphot < 1.7 (likely H/3/[0lll] z ^ 1.4 emitters) and 
2.0 < ^phot < 2.5 (likely z = 2.23 [Oil] emitters) are removed to 
further reduce contamination from higher redshift emitters. Sources 
with spectroscopically confirmed redshifts are included and sources 



" See UKIDSS UDS website for a redshift compilation by O. Almaini and 
the COSMOS data archive for the catalogues, spectra and information on 
the various instruments and spectroscopic programmes. 



with other spectroscopically confirmed lines are removed. In prac- 
tice, 6 Ha spectroscopic sources missed by the selection criteria 
are introduced in the sample; 7 sources found in the sample are 
not Ha - a mix of [Sll], [Nil] and [Olll] emitters. A total of 95 
sources are spectroscopically confirmed as Ha, while 237 sources 
are confirmed as dual Ha-[Olll] emitters. 



2.6.3 Ha emitters at z — 1.47 

Note that, as described in S12, the NBh filter can be combined with 
NB921 (probing the [Oil] emission line), to provide very clean, 
complete surveys of z — 1.47 line emitters, as the filter profiles are 
extremely well-matched for a dual Ha- [Oil] survey. By applying 
the dual narrow-band selection, a total of 346 Ha- [Oil] emitters 
are robustly identified in COSMOS and UDS. However, the dual 
narrow-band selection is only complete if the NB921 survey probes 
down to [Oll]/Ha ^ 0.1 (c.f. SI 2), which is not the case for the 
deepest NBh COSMOS coverage. Additionally, only the central 
1 . 1 deg^ region of the COSMOS field has been targeted with the 
NB921 filter. 

In order to select Ha emitters in areas where the NB921 is not 
deep enough to provide a complete selection, or where NB921 data 
are not available, the following steps are taken. Sources are selected 
if 1.35 < z^hot < 1.55, or if they satisfy the z ^ 1.5 BzK (B — z 
vs. z — K) criteria defined in Figure|4] which is able to recover the 
bulk of the dual narrow-band emitters and sources with high qual- 
ity photometric redshifts of z ^ 1.5. However, the z ^ 1.5 BzK 
selection, although highly complete, is still contaminated by higher 
redshift emitters. In order to exclude likely higher redshift sources 
an additional ziK (i — z vs. z — K\ see SI 2) colour-colour sepa- 
ration is used (see Figure|5]), in combination with rejecting sources 
with Zphot > 1.8. 

The selection leads to a total sample of 5 15 robust Ha emitters 
at z = 1.47, by far the largest sample of Ha emitters at ^ ^ 1.5. 
Comparing the double NB921 and NBh analysis with the colour 
and photo- 2; selection (for sources for which the NB921 data are 
deep enough to detect [Oil]) shows that the colour and photo- 2; se- 
lection by itself results in a contamination of ^ 15 per cent, and 
a completeness of ^ 85 per cent. However, as the double NB921 
and NBh analysis has been used wherever the data are available 
and sufficiently deep, the contamination of the entire sample is es- 
timated to be lower 5 per cent), and the completeness higher 
(^95 per cent). 



2.6.4 Ha emitters at z — 2.23 

As can be seen from the photometric redshift distribution in Fig- 
ure |3] the high quality photo-zs in the COSMOS and UDS fields 
can provide a powerful tool to select z = 2.23 sources. However, 
the sole use of the photometric redshifts can not result in clean, 
high completeness sample of z — 2.23 Ha emitters, not only be- 
cause reliable photometric redshifts are not available for 35 per 
cent of the NBk emitters, at the faint end, but also because the 
errors in the photometric redshifts will be much higher at z ^ 2.2 
than at lower redshift (particularly as one is selecting star-forming 
galaxies). Nevertheless, although spectroscopy only exists for a few 
Ha z = 2.23 sources (zCOSMOS and UDS compilation), double 
line detections between NBk and one of NBh ([Olll]) and/or NB j 
([On]) allow the identification of 84 secure Ha emitters. These can 
be used to optimise the selection criteria and estimate the complete- 
ness and contamination of the sample. 
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Table 3. A summary of the number of sources, narrow-band emitters and Ha-selected emitters for the various surveys undertaken with different narrow-band 
filters. Fields are C: COSMOS, U: UDS. Volumes are presented in units of 10^ Mpc^ and correspond to the total volumes probed; deeper pointings cover 
relatively smaller volumes. SFRs (limiting SFRs) are presented uncorrected for dust extinction. The number of spectroscopically confirmed Ho; emitters are 
presented in the z Ho; conf column, while the number of Ho; sources which are confirmed by an emission-line detection in another narrow-band filter are 
presented in Conf 21ines. Note that the SFR limits represent the limit over the deepest field(s). Also, note that for all four redshifts the number of z -included 
sources missed by the selection is typically < 10, while the 2;-rejected sources are typically < 10 as well. Therefore, decrease of available spectroscopic 
redshifts with redshift does not introduce any bias. 



Filter 


Field 


Detect 


W/ Colours 


Emitters 


Stars 


Spurious 


Ha 


z Ho; conf. 


Conf 21ines 


Volume 


Ha SFR 


NB 


C/U 




# 


3S 


# 


# 


# 


# 


# 


10^ Mpc^ 


Moyr-i 


NB921 


c 


155542 


148702 


2885 (3066) 


247 




521 


38 




5.1 


0.03 


NB921 


u 


236718 


198256 


5799 (6574) 


775 




1221 


8 




5.1 


0.01 


NBj 


c 


32345 


31661 


689 (773) 


40 


46 


425 


81 


158 


7.9 


1.5 


NBj 


u 


21233 


19916 


497 (574) 


49 


30 


212 


14 


79 


11.1 


1.5 


NBh 


c 


65912 


64453 


701 (818) 


60 


63 


327 


28 


158 


68.9 


3.0 


NBh 


u 


26084 


23503 


354 (382) 


23 


5 


188 


18 


188 


22.8 


5.0 


NBk 


c 


65301 


64008 


754 (828) 


54 


23 


346 


4 


65 


31.2 


7.0 


NBk 


u 


28276 


26062 


360 (395) 


32 


10 


172 


2 


19 


22.4 


10.0 


H2 


c 


1054 


940 


52 (57) 


3 


2 


31 





2 


0.9 


3.5 


H2 


u 


1193 


1059 


33 (41) 


7 


1 


14 








0.9 


3.5 



The selection of Ha emitters is done in the same way for both 
COSMOS and UDS, and for both WFCAM and HAWK-I data. An 
initial sample of z — 2.23 Ha emitters is obtained by selecting 
sources for which 1.7 < ^phot < 2.8, where the limits were deter- 
mined using the distribution of photometric redshifts found for con- 
firmed Ha emitters at z = 2.23 (this selects 345 sources, of which 
only 2 are spectroscopically confirmed to be contaminants). Be- 
cause some sources lack reliable photometric redshifts, the colour 
selection [z — K) > (B — z) is used recover additional 2: ~ 2 faint 
emitters. This colour-colour selection is a slight ly modified version 
of the standard BzK colour-colour separation ( |Daddi et al.|2004p" 
It selects 307 additional Ha candidates (and re-selects 331 [96%] 
of those selected through photometric redshifts), and guarantees a 
high completeness of the Ha sample (see Figure |4|. However, the 
BzK selection also selects z ^ 3.3 H/3/[0lll] emitters very ef- 
fectively, and the contamination by such emitters needs to be min- 
imised. In order to do this, sources with 2;phot > 3.0 are excluded 
(64 sources). For sources for which a photometric redshift does not 
exist, a rest-frame UV colour-colour separation is used (B — R vs. 
U—B;see Figure|5] probing the rest-frame UV), capable of broadly 
separating z = 2.23 and z ^ 3.3 emitters due to their different UV 
colours (see Figure|5j this removes a further 29 sources). Four fur- 
ther sources are removed as they are confirmed contaminants (a 
mix of Pa/3, [Slll] and [Olll] at z = 0.65, z = 1.23 and z = 3.23, 
respectively). 

Overall, the selection leads to a total sample of 556 Ha emit- 
ters, by far the largest sample of ^ = 2.23 Ha emitters, and an 
order of magnitude larger than the previous largest samples pre- 
sented by |Geach et al.| (2008 ) and H ayes et all (2010). With the 
limited spectroscopy available, it is difficult to accurately determine 
the completeness and contamination of the sample, but based on the 



^ The selection was modified because the Daddi et al. cut was designed to 
select z > 1.4 sources, while here z = 2.23 emitters are targeted. The 
precise location of the new cut, which is 0.2 magnitudes higher than that of 
Daddi et al., is motivated by the confirmed Ho; emitters, and all of them are 
recovered by the modified cut). 



double-line detections the completeness is estimated to be > 90 per 
cent, and contamination is likely to be 10 per cent. 



3 ANALYSIS AND RESULTS: Ha LF OVER 11 GYRS 

3.1 Removing the contamination by the [Nil] line 

Due to the width of all filters in detecting the Ha line, the ad- 
jacent [Nil] lines can also be detected when the Ha line is de- 
tected at the peak transmission of the filter. A correction for the 
[Nil] line contamination is therefore done, following the rela- 
tion given in SI 2. The relation has been derived to reproduce 
the full SDSS relation between the average log([Nll]/Ha), /, and 
log[EWo([Nll]-hHa)], E: f ^ -0.924 + 4.802^; - 8.892£;2 + 
6.701^^ - 2.27^^ + 0.279^^ This relation is used to correct all 
Ha fluxes at z — 0.4, 0.84, 1.47 and 2.23. The median correction 
(the median [Nll]/([Nll]-hHa)) is ^ 0.25. 

3.2 Completeness corrections: detection and selection 

It is fundamental to understand how complete the samples are as a 
function of line flux. This is done using simulations, as described in 
S09 and further detailed in SI 2. The simulations consider two ma- 
jor components driving the incompleteness: i) the detection com- 
pleteness (which depends on the actual imaging depth and the aper- 
tures used) and ii) the incompleteness resulting from the selection 
(both EW and colour significance). 

The detection completeness is estimated by placing sources 
with a given magnitude at random positions on each individual 
narrow-band image, and studying the recovery rate as a function 
of the magnitude of the source. For the large Subaru frames, 2500 
sources are added for each magnitude, for WFCAM images 500, 
and for HAWK-I frames 100 sources are added for each realisa- 
tion. 

The individual line completeness estimates are performed in 
the same way for the data at the four different redshifts. A set of 
galaxies is defined, which is consistent with being at the approx- 
imate redshift (applying the same photometric redshift -1- colour- 
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Figure 4. Top panels: colour-colour separation of the different emitters in the various samples used to isolate Ho; emitters at the four different redshifts. 
Different symbols show narrow-band emitters with different photometric redshifts. The {B — R) vs. {i — K) colour-colour separation is used to cleanly select 
z = 0.4 Ha emitters with the NB921 data, but also 2; = 0.84 Ho; emitters. The z = OA selection points are [0.1,0.18] and [1.24,2.01], while the selection of 
z = 0.84 Ho; emitters is obtained with the set of points: [—0.26,0.26], [0.1, 0.1 8], [1.0 1,3.21] and [1.24,2.01] (see selection regions shown in the top panels of 
the Figure). Bottom panels: the {B — z) vs {z — K) colours are used in different ways (to maximise completeness) to obtain first pass samples of potential 
z = 1.47 Ha emitters (defined by points: [-0.15,0.40], [0.41,0.0], [2.31,3.6] and [2.41,2.0]) and z = 2.23 Ha emitters jPaddi et al.|2004| z-K> B-z). 
Note, however, that for the Ha samples at z = 1.47 and z = 2.23 the BzK selection is not able to differentiate effectively between Ha and higher redshift 
emitters which can significantly contaminate the samples (see photometric redshift distribution). The higher redshift emitters are filtered out of the samples 
using another set of colours - see Figure|5] 



colour selections to all NB detected sources with no significant ex- 
cess) but not having emission lines above the detection limit. Emis- 
sion lines are then added to the sources, and the study of the recov- 
ery fraction is undertaken. The average completeness corrections as 
a function of Ha flux are presented in Figure |6] Note that the sim- 
ulations include the different EW/colour cuts used in selecting line 
emitters in all bands, and therefore take the EW limits and colour 
selection into account. Also note that because of the very differ- 
ent distributions of magnitudes of Ha emitters from low to higher 
redshift, the EW/colour cut is a much more important source of 
incompleteness for low redshift Ha emitters, and an almost negli- 
gible source of incompleteness at the highest redshift, z — 2.23. 

It should be noted that because of the differences in depth, 
simulations are conducted for each individual frame, and the appro- 



priate completeness corrections applied accordingly when comput- 
ing the luminosity function. For any given completeness correction 
applied, an uncertainty of 20% of the size of the applied correction 
is added in quadrature to the other uncertainties to account for the 
uncertainties in deriving such corrections. 



3.3 Volume 

At 2; = 0.4, the total area surveyed is 1.68 deg^. The NB921 fil- 
ter, centred at 9196 A and with a FWHM of 132 A can probe the 
Ha line (using the top hat approximation) from Zmin = 0.3907 
to Zmax = 0.4108. This means that the narrow-band filter surveys 
an Ha volume of 5.1 x 10^ Mpc^ deg~^. The Ha survey therefore 
probes a total volume of 8.8 x 10^ Mpc^. 
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Figure 5. Top: Colour-colour separation of z = 1.47 Ho; emitters and 
those at higher redshift (z ~ 2.3 H^/[Oiii], z ^ 3.3 [On]); separation is 
obtained by (2;— K) < 5(i—z)— 0.4. Bottom: colour-colour separation of 
z ^3 H/3/[Oiii] emitters from z = 2.23 Ua emitters, given by (B-R) < 
— 0.55(?7— S) + 1.25. Using the BzK colour-colour separation only results 
in some contamination of the Ha sample with higher redshift emitters - 
the B — R \s U — B colour-colour separation allows to greatly reduce 
that. The Figure also indicates the location of each final Ha selected source 
in the colour-colour plot. Note that sources shown as Zphot ~ 1-5 and 
Zphot ~ 2.2 are those with photometric redshifts which are within ±0.2 
of those values. 



The NBj filter (FWHM of 140 A) can be approximated by a 
top hat, probing Zmin = 0.8346 to Zmax = 0.8559 for Ha line 
detections, resulting in surveying 1.5 x 10^ Mpc^ deg~^. As the 
total survey has covered 1.3 deg^, it results in a total volume of 
1.9 X 10^ Mpc^. Assuming the top hat (TH) model for the NBh 
filter (FWHM of 211.1 A, with A™ = 1.606 |^m and A™^ = 
1.627 ^im), the Ha survey probes a (co-moving) volume of 3.3 x 
10^ Mpc^ deg~^. Volumes are computed on a field by field basis 
as each field reaches a different depth (although the difference in 
volume is only important at the faintest fluxes). The total volume 
of the survey is 7.4 x 10^ Mpc"^. The volume down to the deepest 
depth is 3.9 x 10^ Mpc^ (see Table 4 for details). The NBk filter is 
centred on A = 2.121 [im, with a FWHM of 210 A. Using the top 
hat approximation for the filter, it can probe the Ha emission line 
from Zmin 2.2147 to 2;max 2.2467, so with a = 0.016. 
The H2 filter therefore probes a volume of 3.8 x 10^ Mpc^ deg~^. 

The HAWK-I survey uses a slightly different H2 filter, cen- 
tred on A = 2.125 |^m, with FWHM = 300 A. A top hat is an even 
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Figure 6. Average completeness of the various narrow-band surveys as a 
function of Ha flux. Note that the completeness of individual fields/frames 
for each band can vary significantly due to the survey strategy (e.g. see 
difference between one of the deep pointings in COSMOS, NBh D and 
the average over all NBh fields), and thus the completeness corrections are 
computed for each individual fields. 



better approximation of the filter profile, with Zmin = 2.2139 to 
Zma^ = 2.2596 for Ha line detections. The filter effectively probes 
5.5 X 10^ Mpc^ deg"^ Each HAWK-I pointing covers only about 
13.08 arcmin^ and so the complete HAWK-I survey (COSMOS 
and UDS, 0.0156 deg^) probes a total volume of 1.7 x 10^ Mpc^. 
Note that the survey conducted by [Hayes et aL| ( [20T0] ) (using a nar- 
rower NB filter), although deeper, only probed 5.0 x lO^Mpc^, 
so a factor of 3 smaller in volume and over a single field. Table 
4 presents a summary of the volumes probed as a function of Ha 
luminosity and the number of sources detected at each redshift. 



3.4 Filter Profiles: volume corrections 

None of the narrow band filters are perfect top hats (see Figure 
1). In order to model the effect of this bias on estimating the vol- 
ume (luminous emitters will be detectable over larger volumes - 
although, if seen in the filter wings, they will be detected as fainter 
emitters), a series of simulations is done, following S09 and S12. 
Briefly, a top hat volume selection is used to compute a first-pass 
(input) luminosity function and derive the best fit. The fit is used to 
generate a population of simulated Ha emitters (assuming they are 
distributed uniformly across redshift); these are then folded through 
the true filter profile, from which a recovered luminosity function is 
determined. Studying the difference between the input and recov- 
ered luminosity functions shows that the number of bright emit- 
ters is underestimated, while faint emitters can be slightly overesti- 
mated (c.f. S09 for details), but the actual corrections are different 
for each filter and each input luminosity function. This allows cor- 
rection factors to be estimated - these are then used to obtain the 
corrected luminosity function. Corrections are computed for each 
individual narrow-band filter. 



3.5 Extinction Correction 

The Ha emission line is not immune to dust extinction. Measur- 
ing the extinction for each source can in principle be done by sev- 
eral methods, one of which is the comparison between Ha and far- 
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infrared determined SFRs (see Ibar et al. in prep.), while the spec- 
troscopic analysis of Balmer decrements also provides a very good 
estimate of the extinction. As shown in SI 2, the [Oll]/Ha line ratio 
can also be reasonably well calibrated (using Balmer decrement) 
as a dust extinction indicator. For the COSMOS z = 1.47 sample, 
this results in An^ = 0.8 mag (although there is a bias towards 
lower extinction due to the fact that the NB921 survey is not deep 
enough to recover sources with much higher extinctions). However, 
for UDS (where a sufficiently deep NB921 coverage is available) 
a Ahck ~ 1 mag of extinction at Ha is shown to be an appropri- 
ate median correction at z = 1.47 (see S12). That is also similar 
has been found at z = 0.84 ( |Garn et al.|20T0l A^c ~ 1.2). The 
dependence of extinction on observed luminosity is also relatively 
small (S 12) at 2; ^ 1.5 - therefore, for simplicity and for an easier 
comparison, a simple 1 mag of extinction is applied for the four 
redshifts and for all observed luminosities. Note that S12 still find 
a relatively mild luminosity dependence, but one which is offset to 
the local Universe relation (e.g. [Hopkins et al.|2001| ) by 0.5 mag in 
Aua (or, alternatively, the relation is not offset once luminosities 
are divided by at that epoch and compared with local lumi- 
nosities divided by the local ivH^)- 



3.6 Ha Luminosity Functions at z = 0.40, 0.84, 1.47, 2.23 

By taking all Ha selected emitters at the four different redshifts, 
the Ha luminosity function is computed at 4 very different cos- 
mic times, reaching a common observed luminosity limit of ^ 
10^^"^ erg s~^ for the first time in a consistent way over ^11 Gyrs. 
As previously described, the method of S09 and S12 is applied to 
correct for the real profile (see Section [34t . Candidate Ha emitters 
are assumed to be at 2; = 0.4, 0.84, 1.47 and 2.23 for luminosity 
distance calculations. Results can be found in Figure |7] and Table 
|4] Errors are Poissonian, but they include a further 20% of the total 
completeness corrections added in quadrature. 

All derived luminosity functions are fitted with Schechter 
functions defined by three parameters a, 0* and L* : 



0(L)dL = <^* (-^)' 



-(i/i*). 



(t-Y 



(2) 



which are found to provide good fits to the data at all redshifts. In 
the log form, the Schechter function is given by: 



(3) 



0(L)dL = liilO.^* (-^) e-(^/^*> (A)dlogL. 

Schechter functions are fitted to each luminosity function. The best 
fits for the Ha luminosity functions at z = 0.4 — 2.23 are pre- 
sented in Table [5] together with the uncertainties on the parame- 
ters (1 a). Uncertainties are obtained from either the la deviation 
from the best-fit, or the la variance of fits, obtained with a suite of 
multiple luminosity functions with different binning - whichever is 
higher (although they are typically comparable). The best-fit func- 
tions and their errors are also shown in Fi gure [7] to gether with the 
z ^ luminosity function determined bylLy et al.| ( [2007] ) - which 
has extended the work by IGalleg o et al.| {|1995^ at z ^ 0, for a 
local-Universe comparison. Deeper data from the literature are also 
presented for comparison; |Ly et a l. (2011) for z = 0.8, and Hayes 
|et al.| ( [20 10) for 2; = 2.23, after applying the small corrections to 
ensure the extinction corrections are consistent 

The results not only reveal a very clear evolution of the Ha 



^ The correction is applied to obtain data points corrected for extinction by 
1 mag at Ho;. 



luminosity function from z = to z = 2.23, but they also allow 
for a detailed investigation of exactly how the evolution occurs, in 
steps of 2 — 3 Gyrs. The strongest evolutionary feature is the 
increase in L^^ as a function of redshift from z = to z = 2.23 
(see Figure [S]), with the typical Ha luminosity at z ^ 2 (L^^) 
being 10 times higher than locally. This is clearly demonstrated 
in Figure [8] which shows the evolution of the Schechter function 
parameters describing the Ha luminosity function. The evo- 
lution from z ^ to 2; ^ 2.2, can be simply approximated as 
log L* = 0.45^ + log Ll^Q, with log L*^o = 41.87 (see Figure 
[s]). At the very bright end (L > 4L*), and particularly at z > 1, 
there seems to be a deviation from a Schechter function. Follow-up 
spectroscopy of such luminous Ha sources has recently been ob- 
tained for a subset of the z = 1.47 sample, and unveil a significant 
fraction of narrow and broad-line AGN (with strong [Nil] lines as 
well) which become dominant at the highest luminosities (Sobral et 
al. in prep). It is therefore likely that the deviation from a Schechter 
function is being mostly driven by the increase in the AGN activ- 
ity fraction at such luminosities, particularly due to the detection of 
rare broad-line AGN and from very strong [Nil] emission. 

The normalisation of the Ha luminosity function, 0* , is also 
found to evolve, but much more mildly. There is an increase of 0* 
up to 2; ^ 1 (by a factor ^ 4), and then this decreases again for 
higher redshifts by a factor of ^ 2 from z ^ 1 to z = 2.23 - see 
Figure [8] 

The faint-end slope, a, is found to be relatively steep from 
2; up to 2; = 2.23 (when compared to a canonical a = —1.35), 
and it is not found to evolve. The median a over < z < 2.23 is 
-1.60 ±0.08. Very deep data from |Hayes e'tari |20T0^ and |Ly et al.| 
(|2011| not only agree well with such faint-end slope, but even more 
importantly, their data at the faintest luminosities are also very well 
fitted by the best- fit z = 2.23 and z = 0.84 luminosity func- 
tions, even though their data points were not used for the fit. It 
is therefore shown that by measuring the Ha luminosity function 
in a consistent way, and using multiple fields, the faint-end slope 
can be very well approximated by a constant a = —1.6 at least 
up to 2; = 2.23. This shows that while the faint-end slope truly 
is steep at z ^ 2, it does not become significantly steeper from 
2; ^ to 2; ^ 2, and rather has remained relatively constant for the 
last 1 1 Gyrs. The potential steepening of the faint-end slope, which 
has been previously reported (e.g. [Hayes et al.|20"T0| simply results 
from comparing different data- sets which probe different ranges in 
luminosity, use different completeness corrections, different selec- 
tion of emitters and probe a different parameter space. Furthermore, 
the results from Sobral et al. (2011) show that the faint-end slope 
depends relatively strongly on environment, which indicates that 
the changes in the faint-end slope measured before may also have 
resulted by the relatively small areas which can (by chance) probe 
different environments. Note that this is not the case for this pa- 
per because the multi-epoch Ha surveys cover ~ 2 deg^ areas over 
two independent fields and are therefore covering a wide range of 
environments. 

The steep faint-end slope of the Ha luminosity function is in 
very good agreement with the UV luminosity function at 2; 2 and 
above, and particularly consistent with a relatively non-evolving 
a ^ —1.6. This can be seen by comparing the results in this 
paper with those presented by Treye r et al.'('1 998| ), |Arnouts et al.| 
12005 ), and more recently, ||Oesch et al. ( 2010) . It is also likely that 
(similarly to the Ha luminosity function) the large scatter, and the 
different selections/corrections applied, have driven studies to as- 
sume/argue for a steepening of the UV luminosity faint-end slope, 
just like for the Ha luminosity function - see |Oesch et al.| ( [20T0| . 
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Figure 7. The Ho; luminosity function evolution revealed by deep and wide narrow-band surveys at z = 0.4, z = 0.84, z = 1.47 and z = 2.23 presented 
in this work. The data are combined and compared to other studies (Gallego et al. 1995 | Ly et al. 12007) |Geach et al. 12008) [Hayes et al.|2010) , and clearly 
confirm the strong evolution of the luminosity function in the last ~ 11 Gyrs, mostly through a continuous increase in from z = to z = 2.23. Note 
that the z = 0.4 Ho; luminosity function is constrained down to even lower luminosities which are not shown in the figure. Also, note that all Ha luminosities 
have been corrected by extinction with Ahck = 1 mag, and that SFRs shown are based on that extinction correction. SFRs derived directly from observed Ha 
luminosities are a factor of 2.5 lower. 



Overall, the results imply that the bulk of the evolution of the 
star-forming population from z = to z ^ 2.2 is occurring as 
a strong boost in luminosity of all galaxies. The UV luminosity 
results also show very similar trends to the Ha LF, by revealing that 
the strongest evolution to 2; ^ 2 is in the typical luminosity/break 
of the luminosity function, which evolves significantly. However, 
individual measurements for the UV luminosity function at 2; < 
2 are still significantly affected by cosmic variance, small sample 
sizes and much more uncertain dust corrections, and thus the Ha 
analysis provides a much stronger constraint on the evolution of 
star-forming galaxies up to 2; ^ 2.2. 



calibration of Kennicutt (1998) is used to convert the extinction- 
corrected Ha luminosity to a star formation rate: 

SFR(M0yr-') = 7.9 x IQ-^" Lhc (ergs"'), (4) 

which assumes continuous star formation. Case B recombination 
at Te = 10^ K and a Salpeter initial mass function ranging from 
0.1-100 M0. 

As detailed before, a constant 1 magnitude of extinction at Ha 
is assumed for the analysis, which is likely to be a good approach 
for the entire integrated star formation. 



4 THE STAR FORMATION HISTORY OF THE 
UNIVERSE: THE LAST 11 GYRS WITH Ha 

Unveiling the star formation history of the Universe is fundamental 
to understand how, when and at what pace galaxies assembled their 
stellar masses. The best-fit Schechter function fit to the Ha lumi- 
nosity functions at 2; = 0.4, 0.84, 1.47 and 2.23 can be used to esti- 
mate the star formation rate density at the four epochs, correspond- 
ing to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyrs. The standard 



4.1 Removal of the AGN contribution 

Interpreting the integral of the Ha luminosity function as a star 
formation rate density requires a good estimation of the possible 
contribution of AGN to that quantity. For the 2; = 0.84 sample, 
Garn et al.| ( [20T0) ) conducted a detailed search for potential AGN, 
finding a fraction of 8 =b 3% within the Ha population at 2; = 0.84. 
Similar (i.e., ~ 10 %) AGN contaminations at lower redshift have 
also been found by other studies, and therefore assuming a 10% 
contribution from AGN up to 2; ^ 1 is likely to be a good ap- 
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Figure 8. The evolution of the Schechter function parameters which best fit the Ho; luminosity function since z = 2.23. Top left: the evolution of L^^ 
as a function of redshift, revealing that the break of the luminosity function evolves significantly from z = to z = 2.23 by a factor of 10, which can be 
simply parameterised by log L^^ = 0.452; + 41.87. Top right: the evolution of (/>*, which seems to rise mildly up to z ~ 1 and decrease again up to 
z 2.2. Bottom left: the faint end slope, however, is not found to evolve at all from z = 0.0 to z = 2.23 within the scatter and the errors, pointing towards 
a = —1.60 ± 0.08 for the faint-end slope of the Ha luminosity function across the last 11 Gyrs. Bottom right: The 1, 2 and 3 a contours of the best fits to 
the combination of L^^ and fixing a = —1.6. 



Table 5. The luminosity function and star formation rate density evolution for 0.4 < z < 2.2, as seen through a completely self-consistent analysis using 
HiZELS. The measurements are obtained Sit z = 0.4, 0.84, 1.47 and 2.23, correcting for 1 mag extinction at Ho;. Columns present the redshift, break of the 
luminosity function, L"^^, normalisation i4>}j^) and faint-end slope (a) of the Ha luminosity function. The two right columns present the star formation rate 
density at each redshift based on integrating the luminosity function down to ~ 3 Mq yr~i or 41.6 (in log erg s~^) and for a full integration. Star formation 
rate densities include a correction for AGN contamination of 10% at z = 0.4 and z = 0.84 (see |Garn et al. 12010) and 15% at both z = 1.47 and z = 2.23. 
Errors on the faint-end slope a are the 1 a deviation from the best fit, when fitting the 3 parameters simultaneously. As a is very well-constrained at all 
redshifts, and shown not to evolve significantly, L^^ and (/>J^^ are obtained by fixing a = —1.6, and the 1 a errors on L"^^ and cj)^^ are derived from such 
fits (with fixed a). 
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Table 4. Luminosity Functions from HiZELS. Lh^ has been corrected for 
both [Nil] contamination and for dust extinction (using Ah a = 1 mag). 
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proximation. At higher redshifts, and particularly for the sample at 
z = 1.47 and z = 2.23, the AGN activity could in principle be 
different. By looking at a range of AGN indicators - X-rays, radio 
and IRAC colours (and emission lines ratios for sources with such 
informatior[^, it is found that ^ 15 % of the sources are potentially 
AGN at z = 1.47. Similar results are found at z = 2.23. There- 
fore, when converting integrated luminosities to star formation rate 
densities at each epoch, it is assumed that AGNs contribute 10% of 
that up to 2; 1 and 15% above that redshift. While this correction 
may be uncertain, the actual correction will likely be within 5% of 
what is assumed, and in order to guarantee the robustness of the 
measurements, the final measurements include the error introduced 
by the AGN correction - this is done by adding 20% of the AGN 
correction in quadrature to the other errors. 

4.2 The Ha Star Formation History of the Universe 

The results are shown in Table 5, both down to the approximate 
common survey limits, and by fully integrating the luminosity func- 
tion. Figure [9] also presents the results (fully integrating down the 
luminosity functions), and includes a comparison between the con- 
sistent view on the Ha star formation history of the Universe de- 
rived in this paper with the various measurements from the liter- 
ature (Hop kins & Beacom|2006[ |Ly et al.|2011| ), showing a good 
agreement. 

The improvement when compared to other studies is driven 
by: i) the completely self-consistent determinations, ii) the signif- 
icantly larger samples, and iii) the fact that the faint-end slope is 
accurately measured from z ^ to z ^ 2.23 and luminosity func- 
tions determined down to a much lower common luminosity limit 
than ever done before. A comparison with all other previous mea- 
surements (which show a large scatter) reveals a good agreement 
with the Ha measurements. However, the homogeneous Ha anal- 
ysis provides, for the first time, a much clearer and cleaner view of 
the evolution. The results presented in Figure [9] reveal the Ha star 
formation history of the Universe for the last ^11 Gyrs. The evo- 
lution is particularly steep up to about z ^ 1. While the evolution 
is then milder, psFR continues to rise, up to at least z ^ 2. 

Up to z ^ 1, the Ha star formation history is well fitted by 
log psFR = 4x(2; + l) — 2.08. However, such parameterisation 
is not a good fit for higher redshifts. It is possible to fit the entire 
Ha star formation history since z ^ 2.2, or for the last 1 1 Gyrs by 
the simple parameterisation logpsFR = — 0.14T — 0.23, with T 
being the time since the Big Bang in Gyrs (see Figure 9). A power- 
law parameterisation as a function of redshift (ax z)^) yields 
/3 = —1.0, and thus the Ha star formation history can also be 
simply parameterised by log psfr = (ffl) ' clearly revealing that 
PSFK has been declining for the last ~ 11 Gyrs. 

4.3 The Stellar mass assembled in the last 11 Gyrs 

The results presented in this paper can be used to provide an esti- 
mate of the stellar mass (density) which has been assembled by Ha 
star-forming galaxies over the last 1 1 Gyrs. This is done in a similar 
way to |Hopkins & Beacom ( 2006) or |Glazebroo k et al. ( 2004), tak- 
ing into account that a significant part of the mass of newborn stars 

^ Follow up spectroscopy of luminous Ho; sources unveil a significant frac- 
tion of narrow and broad-line AGN which becomes dominant at the highest 
luminosities (Sobral et al. in prep), but is consistent with an overall AGN 
contribution to the Ho; luminosity density of 15 per cent. 
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is indeed recovering essentially all the cosmic star formation happening since z = 2.23. 
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at each redshift is recycled, and can be used in subsequent star for- 
mation episodes. The fraction of recycled mass depends on the IMF 
used. For a Salpeter IMF, which has been used for the Ha calibra- 
tion, the recycling fraction is 30%. Note, however, that changing 
the IMF does not change the qualitative results presented in this 
paper, in particular the agreement between the predicted and the 
measured stellar mass density growth. Nevertheless, changing the 
IMF changes both the normalisation of the star formation history 
and the stellar mass density growth. 

Here, the following approach is taken: the measured stellar 
mass density already in place at 2 ~ 2.2 (many determinations 
exist, e.g. Hopkins ~& Beacoinl|2006[ [Perez-Gonzalez et"aL]|2008[ 
Ilbert et al.|20 09 ) is assumed to be log^o M = 7.45 M© Mpc"^. 
By using the measured Ha star formation history derived in this 
paper (log psFR = — 0.14T — 0.23), a prediction of the evolution 
of the stellar mass density of the Universe is computed, using the 
recycling fraction of the Salpeter IMF (30%). 

The results are presented in Figure[TO| and compared with var- 
ious measurements of the stellar mass density at different redshifts 
available from the literature ([Hopkins & Beacoml[2 006 Perez- 
[Gonzalez et al.|2008 [[Eisner et al.|2008[[Marc hesini'et al. 2009 ). All 
literature results have been converted to a Salpeter IMF if derived 
with a different IMF - including those with a modified Salpeter 
IMF (SalA; resulting in masses a factor of 0.77 lower than Salpeter; 
see e.g. Hopkins & Beacom). 

The results reveal a very good agreement between the predic- 
tions based on the Ha star formation history of the Universe pre- 
sented in this paper since z — 2.23 and the stellar mass density 
evolution of the Universe, measured directly by many authors. The 
results therefore indicate that at least since z — 2.23 the Ha star 
formation history of the Universe is a very good representation of 
the total star formation history of the Universe. It is possible to rec- 
oncile the observed evolution of the stellar mass density with that 
produced from the observed star formation history with very simple 
assumptions, without the need to modify the IMF or have it evolve 
as a function of time. The Ha analysis reveals that star formation 
since z — 2.23 is responsible for 95% of the total stellar mass den- 
sity observed today, with about half of that being assembled from 
z ^ 2.2 io z ^ 1.2, and the other half since z ^ 1.2. Note that 
the same conclusion is reached if the stellar mass density at z = 
is adopted for the normalisation (instead of that at 2; = 2.23), and 
the measured Ha star formation history is used (with appropriate 
recycling factor) to evolve this stellar mass density back to earlier 
epochs. Moreover, if the star formation rate density continues to 
decline with time in the same way as in the last ~ 11 Gyrs, the stel- 
lar mass density growth will become increasingly slower, with the 
stellar mass density of the Universe reaching a maximum which is 
only 5% higher than currently. 



5 CONCLUSIONS 

This paper presents new results from a unique combination of wide 
and deep narrow-band Ha surveys using UKIRT, Subaru and the 
VLT. It has resulted in robust and equally selected samples of sev- 
eral hundreds of Ha emitters in narrow redshift slices, allowing 
to study and parameterise in a completely self-consistent way the 
evolution of the Ha luminosity function over the last 1 1 Gyrs of 
the Universe. The main results are: 

• We robustly select a total of 1742, 637, 515 and 556 Ha emit- 
ters (S>3, EWo(Hc.)>25A) across the COSMOS and the UDS 
fields at 2; = 0.40, 0.84, 1.47 and 2.23, respectively. These are 



by far the largest samples of homogeneously selected Ha emitters, 
while the wide area and the coverage over two independent fields 
allows to greatly overcome cosmic variance and also assemble large 
samples of more luminous galaxies. 

• We find that the Ha luminosity function evolves significantly 
from z^0tO2;^2.2, with the bulk of the evolution being driven 
by the continuous rise in Lhq, by a factor of 10 from the local 
Universe to 2; ^ 2.2, which is well described by log L'^ai^) = 
0.45Z + 41.87 

• By obtaining very deep data over a wide range of epochs, it is 
found that the faint-end slope, a does not evolve with redshift up to 
z - 2.3, and is set to a = -1.60 ± 0.08 for the last 11 Gyrs (0 < 
z < 2.2), contrarily to previous claims (based on heterogeneous 
samples) which argued for a steepening with redshift. 

• The evolution seen in the Ha luminosity function is in good 
agreement with the evolution seen using inhomogeneous compila- 
tions of other tracers of star formation, such as FIR and UV, jointly 
pointing towards the bulk of the evolution in the last 1 1 Gyrs being 
driven by a similar star-forming population across cosmic time, but 
with a strong luminosity increase from z ^ Otoz ^ 2.2. 

• This is the first time Ha has been used to trace SF activity 
with a single homogeneous survey at z = 0.4 — 2.23. The simple 
parametrisations log psfr = — 0.14T — 0.23 (with T being the age 
of the Universe in Gyrs) or log psfr = (z+i) good approxi- 
mations for the last 1 1 Gyrs, showing that psfr has been declining 
since z ^ 2.2. 

• The results reveal that both the shape and normalisation of the 
Ha star formation history are consistent with the measurements of 
the stellar mass density growth, confirming that the Ha cosmic star 
formation history is tracing the bulk of the formation of stars in the 
Universe for ^ < 2.3. 

• The star formation activity over the last ^11 Gyrs is responsi- 
ble for producing ^ 95% of the total stellar mass density observed 
locally today, with about half of that being assembled from ^ 2.2 
to 2; 1.2, and the other half at 2; < 1.2. 

The results presented in this paper provide a self-consistent 
view that improves our understanding of the evolution of star- 
forming galaxies. Particularly, it shows that the evolution of the 
star-forming population in the last 1 1 Gyrs has been mostly driven 
by a change in the typical star formation rate of the population 
(Lhc), while the faint-end slope of the Ha LF has remained con- 
stant, and the change in the normalisation has been much more 
moderate. The strong evolution in L^^ (or SFR*) may well be 
unveiling something very fundamental about the evolution of star- 
forming galaxies, as it seems to mark a transition between disk 
and mergers (e.g. S09) in the last 9-10 Gyrs. Also, scaling SFRs 
by the SFR* at each epoch (or Ha luminosities by L^^ at each 
epoch) seems to recover relatively non-evolving relations between 
scaling SFRs/luminosities and e.g. dust extinction (S12), morpho- 
logical class (S09), or the typical dark matter halo in which the 
star-forming galaxies are likely to reside f Sobral et al.|201 0). 

The results presented in this paper also complement the cur- 
rent view on the evolution of the stellar mass function over the 
last 11 Gyrs (e.g. [Ilbert et al.| [20091 [Peng et al.[|20T0l [March-| 
[esini et al.|2012| ), which also reveal a non-evolving faint-end slope 
(of the stellar mass function) at least for 2; < 2, but shallower, 
a = —1.3. However, the typical mass of the stellar mass func- 
tion, M*, is found to be roughly constant in the last 11 Gyrs, with 
the main change being 0*, which continuously increases in the last 
^ 11 Gyrs. Combining the results of the evolution of the Ha lumi- 
nosity function with those of the evolution of the stellar mass func- 
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tion point towards a star-forming population mostly evolving with 
an overall decrease in their SFRs/luminosity, and with the over- 
all population presenting a typical non-evolving mass and simply 
evolving by a change in number density, consistent with the sce- 
nario proposed by |Peng et al.|p010| . 



ACKNOWLEDGMENTS 

DS is supported by a NOVA fellowship. IRS acknowledges a Lev- 
erhulme Senior Fellowship. PNB acknowledges support from the 
Leverhulme Trust. YM, JPS and IRS thank the U.K. Science and 
Technology Facility Council (STFC). JEG is supported by a Bant- 
ing Fellowship, administered by the Natural Science and Engineer- 
ing Research Council of Canada. We would like to thank Richard 
Ellis, Simon Lilly, Peter Capak, Adam Muzzin, Taddy Kodama, 
Masao Hayashi, Andy Lawrence, Joop Schaye, Marijn Franx, Huub 
Rottgering, Ry chard Bouwens and Renske Smit for many interest- 
ing and helpful discussions. We would also like to thank Chun Ly 
for helpful comments. The authors wish to thank all the JAC staff 
for their help conducting the observations at the UKIRT telescope, 
and their continuous support and we are also extremely grateful 
to all the Subaru staff. We also acknowledge ESO and Subaru 
for service observations. Finally, the authors fully acknowledge 
the tremendous work that has been done by both COSMOS and 
UKIDSS UDS/SXDF teams in assembling such large, state-of-the- 
art multi-wavelength data-sets over such wide areas, as those have 
been crucial for the results presented in this paper. 



REFERENCES 

Arnouts S., Schiminovich D., Ilbert O., Tresse L., Milliard B., 
Treyer M., BardeUi S., et al. 2005, ApJL, 619, L43 

Atek H., Malkan M., McCarthy R, Teplitz H. I., Scarlata C, Siana 
B., Henry A., et al. 2010, ApJ, 723, 104 

Bertin E., Arnouts S., 1996, A&AS, 117, 393 

Best P, Small I., Sobral D., Geach J., Garn T, Ivison R., Kurk J., 
Dalton G., Cirasuolo M., Casah M., 2010, arXiv: 1003.5 183 

Bunker A. J., Warren S. J., Hewett P C, Clements D. L., 1995, 
MNRAS,273, 513 

Casali M., Pirard J.-F., Kissler-Patig M., Moorwood A., Bedin L., 
Biereichel P., Delabre B., et al. 2006, in Society of Photo- Optical 
Instrumentation Engineers (SPIE) Conference Series Vol. 6269 
of Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, HAWK-I: the new wide-field IR imager for 
the VLT 

Cirasuolo M., McLure R. J., Dunlop J. S., Almaini O., Foucaud 
S., Simpson C, 2010, MNRAS, 401, 1166 

Cutri R. M., Skrutskie M. F, van Dyk S., Beichman C. A., Car- 
penter J. M., Chester T, Cambresy L., et al. 2003, 2MASS All 
Sky Catalog of point sources. 

Daddi E., Cimatti A., Renzini A., Fontana A., Mignoh M., 
Pozzetti L., Tozzi P, Zamorani G., 2004, ApJ, 617, 746 

Dutton A. A., van den Bosch F. C, Dekel A., 2010, MNRAS, 405, 
1690 

Eisner F, Feulner G., Hopp U., 2008, A&A, 477, 503 
Gallego J., Zamorano J., Aragon- Salamanca A., Rego M., 1995, 
ApJL, 455, LI 

Garn T, Sobral D., Best P N., Geach J. E., Small I., Cirasuolo 
M., Dalton G. B., Dunlop J. S., McLure R. J., Farrah D., 2010, 
MNRAS, 402, 2017 



Geach J. E., Simpson C, Rawlings S., Read A. M., Watson M., 

2008, VizieR Online Data Catalog, 838, 11369 
Geach J. E., Small I., Best P N., Kurk J., Casah M., Ivison R. J., 

Coppin K., 2008, MNRAS, 388, 1473 
Gilbank D. G., Baldry I. K., Balogh M. L., Glazebrook K., Bower 

R. G., 2010, MNRAS, 405, 2594 
Glazebrook K., Abraham R. G., McCarthy P. J., Savagho S., Chen 

H.-W, Crampton D., Murowinski R., et al. 2004, Nature, 430, 

181 

Gonzalez V, Labbe I., Bouwens R. J., Illing worth G., Franx M., 

KriekM., Brammer G. B., 2010, ApJ, 713, 115 
Hayes M., Schaerer D., Ostlin G., 2010, A&A, 509, L5-h 
Hopkins A. M., 2004, ApJ, 615, 209 
Hopkins A. M., Beacom J. F, 2006, ApJ, 651, 142 
Hopkins A. M., Connolly A. J., Haarsma D. B., Cram L. E., 2001, 

AJ, 122, 288 

Hopkins A. M., Connolly A. J., Szalay A. S., 2000, AJ, 120, 2843 
Ilbert O., Capak P., Salvato M., Aussel H., McCracken H. J., 

Sanders D. B., Scoville N., et al. 2009, ApJ, 690, 1236 
Lawrence A., Warren S. J., Almaini O., Edge A. C, Hambly N. C, 

Jameson R. F, Lucas P, et al. 2007, MNRAS, 379, 1599 
Lilly S. J., Le Brun V, Maier C, Mainieri V, Mignoli M., Scodeg- 

gio M., Zamorani G., et al. 2009, ApJS, 184, 218 
Lilly S. J., Le Fevre O., Hammer F, Crampton D., 1996, ApJL, 

460, LI 

Ly C, Lee J. C, Dale D. A., Momcheva I., Salim S., Staudaher 

S., Moore C. A., Finn R., 2011, ApJ, 726, 109 
Ly C, Malkan M. A., Kashikawa N., Shimasaku K., Doi M., Na- 

gao T, lye M., et al. 2007, ApJ, 657, 738 
Malkan M., Teplitz H., McLean I., 1995, ApJL, 448, L5 
Marchesini D., Stefanon M., Brammer G. B., Whitaker K. E., 

2012, arXiv: 1201.6365 
Marchesini D., van Dokkum P. G., Forster Schreiber N. M., Franx 

M., Labbe I., Wuyts S., 2009, ApJ, 701, 1765 
McCarthy P J., Yan L., Freudling W, Teplitz H. I., Malumuth 

E. M., Weymann R. J., Malkan M. A., et al. 1999, ApJ, 520, 548 
Moorwood A. F. M., van der Werf P P, Cuby J. G., Oliva E., 2000, 

A&A, 362, 9 

Oesch P A., Bouwens R. J., Carollo C. M., Illingworth G. D., 
Magee D., Trenti M., Stiavelli M., Franx M., Labbe I., van 
Dokkum P G., 2010, ApJL, 725, L150 

Ono Y., Ouchi M., Shimasaku K., Akiyama M., Dunlop J., Farrah 
D., et al.. 2010, MNRAS, 402, 1580 

Ouchi M., Shimasaku K., Akiyama M., Simpson C, Saito T, 
Ueda Y, Furusawa H., et al. 2008, ApJS, 176, 301 

Ouchi M., Shimasaku K., Furusawa H., Saito T, Yoshida M., 
Akiyama M., Ono Y, et al. 2010, ApJ, 723, 869 

Ouchi M., Shimasaku K., Okamura S., Furusawa H., Kashikawa 
N., Ota K., Doi M., et al. 2004, ApJ, 611, 660 

Peng Y.-j., Lilly S. J., Kovac K., Bolzonella M., Pozzetti L., Ren- 
zini A., Zamorani G., et al. 2010, ApJ, 721, 193 

Perez-Gonzalez P G., Rieke G. H., Villar V, Barro G., Blaylock 
M., Egami E., Gallego J., et al. 2008, ApJ, 675, 234 

Pirard J.-F, Kissler-Patig M., Moorwood A., Biereichel P., De- 
labre B., Dorn R., Finger G., et al. 2004, in A. F M. Moor- 
wood & M. lye ed., Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series Vol. 5492 of Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Se- 
ries, HAWK-I: A new wide-field 1- to 2.5-/im imager for the 
VLT. pp 1763-1772 

Salim S., Rich R. M., Chariot S., Brinchmann J., Johnson B. D., 
Schiminovich D., Seibert M., et al. 2007, ApJS, 173, 267 



18 D. Sobral et al. 



Scoville N., Aussel H., Brusa M., Capak R, Carollo C. M., Elvis 

M., Giavalisco M., et al 2007, ApJS, 172, 1 
Shim H., Chary R.-R., Dickinson M., Lin L., Spinrad H., Stern 

D., Yan C.-H., 2011, arXiv: 1103.4124 
Shioya Y, Taniguchi Y, Sasaki S. S., Nagao T., Murayama T., 

Takahashi M. I., Ajiki M., et al. 2008, ApJS, 175, 128 
Simpson C, Martinez-Sansigre A., Rawlings S., Ivison R., 

Akiyama M., Sekiguchi K., Takata T., Ueda Y, Watson M., 

2006, MNRAS, 372, 741 

Small L, Sharp R., Swinbank A. M., Akiyama M., Ueda Y, Fou- 

caud S., Almaini O., Croom S., 2008, MNRAS, 389, 407 
Sobral D., Best R, Matsuda Y, Small I., Geach J., Cirasuolo M., 

2012,arXiv:1109.1830 
Sobral D., Best R N., Geach J. E., Small L, Cirasuolo M., Garn 

T, Dalton G. B., Kurk J., 2010, MNRAS, 404, 1551 
Sobral D., Best R N., Geach J. E., Small I., Kurk J., Cirasuolo M., 

Casah M., et al. 2009, MNRAS, 398, 75 
Sobral D., Best R N., Small L, Geach J. E., Cirasuolo M., Garn 

T., Dalton G. B., 2011, MNRAS, 411, 675 
Straughn A. N., Kuntschner H,, Kummel M., Walsh J. R., Cohen 

S. H., Gardner J. R, Windhorst R. A., et al. 2011, AJ, 141, 14 
Tadaki K.-L, Kodama T., Koyama Y, Hayashi M., Tanaka I., 

Tokoku C, 201 1, PASJ, 63, 437 
Treyer M. A., Elhs R. S., Milhard B., Donas J., Bridges T. J., 

1998, MNRAS, 300, 303 
van Breukelen C, Cotter G., Rawlings S., Readhead T., Bonfield 

D. , Clewley L., Ivison R., Jarvis M., Simpson C, Watson M., 

2007, MNRAS, 382, 971 

van Dokkum P. G., Brammer G., Fumagalli M., Nelson E., Franx 

M., Rix H.-W, Kriek M., Patel S., et al. 2011, ApJL, 743, L15 
Villar v., Gallego J., Perez-Gonzalez R G., Pascual S., Noeske K., 

Koo D. C, Barro G., Zamorano J., 2008, ApJ, 677, 169 
Yagi M., Kashikawa N., Sekiguchi M., Doi M., Yasuda N., Shi- 

masaku K., Okamura S., 2002, AJ, 123, 66 
Yamada T., Kodama T., Akiyama M., Furusawa H., Iwata I., Ka- 

jisawa M., lye M., Ouchi M., et al. 2005, ApJ, 634, 861 
Yan L., McCarthy R J., Freudhng W, Tephtz H. I., Malumuth 

E. M., Weymann R. J., Malkan M. A., 1999, ApJL, 519, L47 



APPENDIX A: CATALOGUES OF CANDIDATE HIZELS 
NARROW-BAND EMITTERS 

The catalogues of potential narrow-band emitters over the COS- 
MOS and UDS fields are presented in Tables A.l (NB921), A.2 
(NBj), A.3 (NBh) and A.4 (NBk). It contains IDs (including 
field and observing band). Right Ascension (RA), Declination 
(Dec), narrow-band magnitude (NB), broad-band magnitude (BB), 
the significance of the narrow-band excess (S), Estimated Flux 
(log 10), estimated observed EW, and a flag for those that are clas- 
sified as Ha. Note that only the online version contains the full cat- 
alogue - here only five entries of the table are shown as examples 
of the entire catalogues. 

This paper has been typeset from a T^C/ file prepared by the 
author. 
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Table Al. A catalogue of all E > 3 NB921 narrow-band sources selected in the UDS and COSMOS fields from HiZELS. NB magnitudes are NB921 (AB); 
BB are z' magnitudes (AB). Note that on the printed version this only contains 5 entries, in order to provide some examples contained in the full catalogue. 
The full catalogue is available on-line. 



ID 


R.A. 


Dec. 


NB 


BB 


S 


log Flux 


EWobs 


Ha selected 




(J2000) 


(J2000) 


(AB) 


(AB) 




ergs"-*^ 


A 




HiZELS-COSMOS-NB921-S12-14 


095813.57 


+02 17 15.6 


23.86+0.11 


24.85+0.08 


13.6 


-16.255 


418.4 


No 


HiZELS-COSMOS-NB921-S12-105 


095957.96 


+021741.5 


23.14+0.03 


23.39+0.02 


9.9 


-16.132 


119.8 


No 


HiZELS-COSMOS-NB921-S12-136 


0958 57.55 


+021745.0 


21.64+0.01 


22.11+0.01 


56.5 


-15.636 


105.9 


Yes 


HiZELS-COSMOS-NB921-S12-4064 


095920.49 


+022054.5 


22.06+0.01 


22.21+0.01 


14.3 


-16.084 


39.13 


No 


HiZELS-UDS-NB921-S 12-220582 


021619.37 


-0513 54.2 


22.37+0.01 


22.65+0.01 


55.6 


-16.287 


38.22 


Yes 



Table A2. A catalogue of all S > 3 NBj narrow-band sources selected in the UDS and COSMOS fields from HiZELS. NB magnitudes are NBj (Vega); 
BB are J magnitudes (Vega). Note that on the printed version this only contains 5 example entries, in order to provide some examples contained in the full 
catalogue. The full catalogue is available on-line. 



ID 


R.A. 


Dec. 


NB 


BB 


S 


log Flux 


EWobs 


Ha selected 




(J2000) 


(J2000) 


(Vega) 


(Vega) 




ergs"-*- 


A 




HiZELS-COSMOS-NBJ-S 12-206 


095842.11 


+015840.5 


21.58+0.20 


22.83+0.23 


3.5 


-16.11 


463.5 


No 


HiZELS-COSMOS-NBJ-S 12-292 


095842.65 


+0227 24.2 


20.51+0.08 


20.85+0.05 


3.7 


-16.08 


65.4 


Yes 


HiZELS-COSMOS-NBJ-S 12-293 


095842.63 


+021953.2 


20.84+0.10 


21.46+0.07 


4.4 


-16.01 


139.3 


Yes 


HiZELS-UDS-NB J-S 12-5 


021616.81 


-05 0909.2 


20.51+0.08 


20.85+0.05 


5.4 


-15.89 


194.6 


Yes 


HiZELS-UDS-NBJ-S 12-26 


021617.15 


-05 0742.4 


20.84+0.10 


21.46+0.07 


4.6 


-15.96 


79.8 


Yes 



Table A3. A catalogue of all S > 3 NBh narrow-band sources selected in the UDS and COSMOS fields from HiZELS. NB magnitudes are NBh (Vega); 
BB are H magnitudes (Vega). Note that on the printed version this only contains 5 example entries, in order to provide some examples contained in the full 
catalogue. The full catalogue is available on-line. 



ID 


R.A. 


Dec. 


NB 


BB 


S 


log Flux 


EWobs 


Ha selected 




(J2000) 


(J2000) 


(Vega) 


(Vega) 




ergs"-*^ 


A 




HiZELS-COSMOS-NBH-S 12-206 


095747.31 


+020208.1 


21.58+0.14 


20.99+0.05 


3.1 


-15.98 


251.4 


No 


HiZELS-COSMOS-NBH-S 12-292 


095747.33 


+015153.1 


20.51+0.02 


18.33+0.01 


12.1 


-15.40 


77.3 


No 


HiZELS-COSMOS-NBH-S 12-293 


095747.65 


+02 16 19.8 


20.84+0.08 


19.80+0.02 


3.3 


-15.87 


101.8 


No 


HiZELS-UDS-NBH-S 12-4012 


021644.17 


-044453.0 


20.04+0.18 


20.93+0.04 


3.0 


-15.87 


308 


Yes 


HiZELS-UDS-NBH-S12-4136 


021645.35 


-0407 26.8 


19.15+0.08 


19.76+0.02 


5.2 


-15.64 


171 


Yes 



Table A4. A catalogue of all S > 3 NBk narrow-band sources selected in the UDS and COSMOS fields from HiZELS. NB magnitudes are NBk (Vega); 
BB are K magnitudes (Vega). Note that on the printed version this only contains 5 example entries, in order to provide some examples contained in the full 
catalogue. The full catalogue is available on-line. 



ID 


R.A. 


Dec. 


NB 


BB 


E 


log Flux 


EWobs 


Ha selected 




(J2000) 


(J2000) 


(Vega) 


(Vega) 




ergs"-*- 


A 




HiZELS-COSMOS-NBK-S 12-2227 


09 58 14.05 


+023405.1 


20.26+0.23 


21.68+0.16 


3.1 


-16.31 


728 


Yes 


HiZELS-COSMOS-NBK-S 12-45390 


100059.58 


+024435.9 


19.24+0.09 


19.92+0.04 


5.4 


-16.09 


208 


No 


HiZELS-COSMOS-NBKS 1 2-44 1 93 


100055.39 


+015955.1 


19.02+0.06 


19.68+0.02 


7.6 


-16.01 


198 


Yes 


HiZELS-UDS-NBK-S12-15961 


021809.03 


-044749.9 


17.47+0.02 


18.27+0.01 


30.6 


-15.33 


262 


Yes 


HiZELS-UDS-NBK-S12-22618 


0218 53.02 


-05 0107.7 


19.36+0.10 


20.02+0.02 


4.5 


-16.16 


196 


No 



